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Abstract 

We present several recent analyses of Dalitz plots from the CLEO-c ex- 
periment, including published and preliminary analyses of D + — > h~tt + tt + , 
D + — > K~ir + ir + , and D° — ► Kg l tt + tt~ decays. More information on these 

analyses can be found in References HI 121 131) New preliminary analyses we 
present include a search for CP asymmetry in D + — > K + K~ir + decays and a 
Dalitz plot analysis of D° -> K%ir°ir°. 

We report on a search for the CP asymmetry in the singly Cabibbo- 
suppressed decay D + — > K + K~ir + using a data sample of 572 pb _1 accumu- 
lated with the CLEO-c detector and taken at the e+e~ — * V(3770) resonance. 
We have searched for CP asymmetries using a Dalitz plot based analysis that 
determines the amplitudes and relative phases of the intermediate states. 

We also use a 281 pb _1 CLEO-c data sample taken at the e + e~ — > 
V>(3770) resonance to study the D° ->■ K^w tt Dalitz plot. Our nominal fit 
includes the K°, K*(892), / (980), / (137Q), and A"* (1680) resonances. 



1 Search for CP asymmetry in D + — > K + K ir + Decays 

Singly Cabibbo-suppressed (SCS) D-meson decays are predicted in the Stan- 
dard Model (SM) to exhibit CP-violating charge asymmetries smaller than the 
order of 10 -3 . Direct CP violation in SCS decays could arise from the interfer- 
ence between tree-level and penguin processes. Doubly Cabibbo-suppressed 
and Cabibbo-favored (CF) decays are expected to be CP invariant in the 
SM due to the lack of contribution from penguin processes. Measurements 
of CP asymmetries in SCS processes greater than 0(1O~ 3 ) would be evidence 
of physics beyond the SM © , 

We define two variables: the energy difference AE = J^. Ei — E^ eam and 

the beam-constrained mass mec = "J ^ cam — | J2% Pi\ 2 i where Ei, Pi are the 
energy and momentum of each D decay product, and -Ebcam is the beam energy. 
We define a signal box corresponding to 2.5 standard deviations in each vari- 
able, and remove multiple candidates in each event by choosing the candidate 
that gives the smallest \AE\. We obtain 13693 ± 137 D+ K+R-TT+ sig- 
nal candidates. To reduce smearing effects introduced by the detector, a mass 
constraint fit for the D + candidate is applied to obtain the mass squared vari- 
ables, m 2 K+ + and m 2 K _ 7T+ , for the D + — > K + K~ir + Dalitz plot (DP) shown 
in Figure [T|a). 

The decay amplitude as a function of DP variables is expressed as a 
sum of two-body matrix elements and one non-resonant (NR) decay amplitude 
For most resonances, the matrix element is parameterized by Breit-Wigner 
shapes that take into account D meson and intermediate resonance form factors 
and angular dependence. For the /o(980) we use a Flatte function E*. For 
the ao(980), we use the function in Ref. We choose the same phase 

conventions for the intermediate resonances as the E687 Collaboration 0) A 
fit fraction (FF), the integral of a single component divided by the sum of 
all components, is reported for each intermediate resonance to allow for more 
meaningful comparisons between results. 

For D + decays to K~n + S- wave states, we consider three amplitude 
models. One model uses a coherent sum of a uniform non-resonant term and 
Breit-Wigner term for the ^(1430) resonance. The second model only uses 
a Breit-Wigner term for the -RTq(1430) resonance. The third model uses the 
LASS amplitude for K~it + — > K~tt + clastic scattering 1*3 GUI). We present 



results only for the third model, although the first model provides a similar fit. 




Figure 1: The results of fitting the D + — > K + K~ir + data for Model three, 
(a) The scatter plot for squared mass of K~tt + versus K + ir + and the projec- 
tions onto squared mass of (b) K~tt + (c) K + tt + and (d) K + K~ for both fit 
(curve) and data (points) are shown. The dashed line shows the background 
contribution. 

We determine the detection efficiency as a function of the two DP variables 
by fitting a signal MC sample generated with a flat distribution in the phase 
space. We use a fit to the events in the AE sideband (24 < \AE\ < 42 MeV 
and \m-QC — m D+\ < 9 MeV/c 2 ) to describe the background distribution of 
the DP. Having information for both the background and efficiency, as well as 
the fraction of signal events in the signal region, we fit the data in the DP to 
extract the amplitudes and phases of any contributing intermediate resonances. 
We perform an unbinned maximum likelihood fit. The signal fraction / is 
fa = (84.1 ± 0.2)%, constrained in the fit to be within its error af obtained 
from the fit to the m-Qc distribution. We begin by fitting the DP with all 
known resonances that may possibly contribute to this decay. We determine 



Table 1: The fit results in Model three. The errors shown are statistical, 
experimental systematic, and modeling systematic respectively. 



Component 



Amplitude 



Phase (°) 



Fit Fraction (%) 



K (892)°Jf+ 
K-ir+(S)K+ 
a (980)7r+ 
0(lO2O)7r + 
/ 2 (1270)tt+ 
a (1450)7r+ 
0(168O)tt+ 



1 (fixed) 



4.53±0.16^ni 3 

0. 74±0.09l^l°;3 6 9 
i oq _i_ n n9+ 00 + 01 

1. zo zsz u.u^_ oo-0 .02 

01 ± 1O+0.03+0.11 

u.yi in u.j.o_q .01-0.24 

1 Q 6±0 iQ+0.20+0.45 



(fixed) 
21 ± 3±°±£ 
96 ±7+^ 

— 148 ± 3±i±l 
20 ± 6±g±?i 

116 ±5™ 



2.6 ± 0.3^;™ -96 ± 10i?+^ 2 



Q c _|_ 1 n+1.6+1.6 
0.0 in i-U_o.o-2.6 



-156 ± 6±S±|° 



23.9±0.6±°'™ 
53 ± 3±g±§ 

174-0 A+0-1 + 1.3 
1. 1 3= U.*_o. _o.6 

28.0 ± 0.5±g:2 ± 0.5 
0.9 ± 0.2±° j ± 0.2 



3.4±0.5±^; 5 2 
i.89±0.18™° ; ; 

9 1-1-1 0+2.4+2.2 
z.i zsz J--^_o.O-1.3 



which resonances are to be included by maximizing the fit confidence level 
(C.L.). The procedure is to add all possible resonances, then subsequently 
remove those which do not contribute significantly, or worsen our C.L. The 
projections of the DP for the fit to Model three are shown in Figures [ljb-d) . 
The results of the fit amplitudes, phases, and fractions including errors are 
shown in Table Q] for Model three. 



Table 2: Acp for each component of the fit using D ± samples in Model three. 
The errors for fit fractions and phases are statistical only, and those for Aqp 
are statistical, experimental systematic, and modeling systematic respectively. 



Component j 



K (892)°i^ 

a (980)7r 4 
0(1O2O)tH 
/ 2 ( 1270)tH 
a (1450)7r 4 
0(1680)7^ 
#2(1430)°^ 



.4 



CP] 



(%) 



-0.1 ±2.9™* 

-i±5+ii 6 4 

-ll±23±4+6 4 
o n + i q+0.1+0.2 

O.U zsz J--»_o.2-0.3 

4 ± 25±%t% 
-18±14t^ 6 
-9 ± 21+4_+ 7 



69 ±51 



+1+8 

-28-41 



To search for CP violation in this model, we fit the D + and D~ sam- 
ples independently. We use the same background fraction and PDF as those 
used in the fit to the total sample, but different coefficients for efficiency func- 
tions which are obtained from signal MC of D decays. The calculated CP 
asymmetry, Acpj 



pp. 

— r D + 



is shown for each resonance j in Table |ij 



2 Dalitz Plot Analysis of D° -> A^7r°7r° Decays 

The PDG nas little information on the D° — > X^7r 7r° decay. In addition 
to providing a more comprehensive study of the D° — > A!j7r 7r decay, this DP 
analysis seems like a good place to look for the low mass irir S-wave signa- 
ture of the a. The Aij7r + 7r~ mode is much cleaner and has better statistics, 
but the p° resonance overlaps the region where we would expect to find the 
low mass S'-wave signature. Using CLEO-c data, we eliminate nearly all of 
the background by doing a double-tagged analysis, where both D mesons are 
completely reconstructed. 

We have analyzed 281 pb _1 of CLEO-c data taken on the e + e~ — > ?/>(3770) 
resonance. In a double-tagged analysis, both D mesons are reconstructed. For 
our double-tagged analysis, we consider candidates with one D reconstructed as 
Kgir°ir°, and the other D reconstructed using any of the following decay modes 
(charge conjugation is implied throughout this analysis): D° — > K + ir~, D° — > 
K + Tc~Tr°, D° — » K + 7r~7r + 7r~. In a single-tagged analysis, we reconstruct only 
one D meson in the event, which decays to KgTr°ir°. 

Table 3: D° — » i^^7r°7r° signal yield, number of candidates, and signal fraction 



Result 


Double Tag 


Single Tag 


Signal Yield 
Total Candidates 
Signal Fraction 


257 ± 17 
276 

0.931 ± 0.062 


1884 ± 56 
2548 
0.739 ± 0.022 



To reduce 2tt background that fakes a Kg, we enforce a 2a enhanced 
flight significance selection criteria on our Kg candidates. To reduce the Ktttt 
background, we require \dE / dx P i on \ < 3a and dE/dxkaon < —2a for both A'lj 
daughter pions. We use the same particle identification selection criteria for 
double-tagged and single-tagged analyses. We apply a 2a selection criteria on 
the reconstructed Kg mass. After enforcing our selection criteria on the Kg 
mass, we apply a 2a selection criteria on AE. We additionally apply a 2a 
cut on the beam constrained mass. For each event that has more than one 
candidate, we require the following: For the double-tagged data, we take the 
average of the signal beam constrained mass and the tagged beam constrained 
mass, and we select whichever candidate's average is closest to the nominal D 



mass. For the single-tagged data, we select the candidate with AE closest to 
zero. Table [3] shows our signal yield and signal fraction. 

For this analysis, we define our DP variables as follows: x = larg er in ts-q , 

s 

y = to^o^o, z = smaller m 2 K0Tr0 . When fitting such a Dalitz plot, we must take 
into account the fact that the two ir° final state particles are indistinguishable, 
so we explicitly symmetrize the functions we use in x and z. 

To study the efficiency of reconstructing our signal, we generate 100000 
signal Monte Carlo events distributed uniformly across the Dalitz plot phase 
space. Half of these events force the D° to decay directly into K s ir°ir° and 
the D to decay into neutrinos. The other half of these events force the D 
to decay directly into our signal mode and the D° to decay into neutrinos. 
We fit the efficiency over the Dalitz plot to a third-order polynomial explicitly 
symmetric in x and z. To fit for the background, we use a sideband from single- 
tagged data which is centered 5<r mD0 lower in tubc than the signal region, with 
the same width as that of the signal region, and has the appropriate range in 
AE which conserves the boundaries of the signal DP. We use this background 
shape for the double-tagged data as well as for the single-tagged data. We fit 
the background events to a third-order polynomial explicitly symmetric in x 
and z, plus a non-interfering if* (892) Breit-Wigner in both x and z. 

The signal is parameterized with an isobar model that has four interfer- 
ing resonances plus one non-interfering resonance. To enforce the symmetry 
requirement in the DP, we include each K* resonance as an x resonance and a 
z resonance, while using the same amplitude and phase for the x contribution 
and z contribution. The parameters for the K s , if* (892), and K* (1680) come 
from the PDG . The parameters for the Jo (980) are approximated from a 
BES paper The parameters for the /q(1370) come from Reference . 

Figure [2ja) displays the DP from the double-tagged data. To fit this 
DP with an unbinned maximum likelihood fitter, we fix the signal fraction to 
0.931 as determined from the beam constrained mass distribution. The fit also 
fixes the efficiency parameters and background parameters as determined from 
the signal Monte Carlo and sideband. The fit determines the amplitudes and 
phases of the resonances and calculates the fit fractions. Figure [2jt>) shows the 
fit results. 

To estimate systematic errors, we use the technique developed by Jim 
Wiss and Rob Gardner . Using this technique, the systematic errors are 




Figure 2: (a) Dalitz plot of the double-tagged data (2 entries per candidate, sec- 
ond entry has x and z swapped) and (b) fits to the double-tagged x+z projection 
and double-tagged y projection. 



essentially independent of the number of systematic sources considered ^LM . 
Table [4] gives our preliminary results. We are currently extending our analysis 
to the full available CLEO-c ^(3770) data sample, and studying the effects of 
using a o or k S- wave to possibly improve our fit. 
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